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ABSTRACT 

Context. Stellar dynamo processes can be explored by measuring the magnetic field. This is usually obtained using the atomic 
and molecular Zeeman effect in spectral lines. While the atomic Zeeman effect can only access warmer regions, the use of 
molecular lines is of advantage for studying cool objects. The molecules MgH, TiO, CaH, and FeH are suited to probe stellar 
magnetic fields, each one for a different range of spectral types, by considering the signal that is obtained from modeling various 
spectral types. 

Aims. We have analyzed the usefulness of different molecules (MgH, TiO, CaH, and FeH) as diagnostic tools for studying 
stellar magnetism on active G-K-M dwarfs. We investigate the temperature range in which the selected molecules can serve as 
indicators for magnetic fields on highly active cool stars and present synthetic Stokes profiles for the modeled spectral type. 
Methods. We modeled a star with a spot size of 10% of the stellar disk and a spot comprising either only longitudinal or only 
transverse magnetic fields and estimated the strengths of the polarization Stokes V and Q signals for the molecules MgH, TiO, 
CaH, and FeH. We combined various photosphere and spot models according to realistic scenarios. 

Results. In G dwarfs, the molecules MgH and FeH show overall the strongest Stokes V and Q signals from the starspot, whereas 
FeH has a stronger Stokes V signal in all G dwarfs, with a spot temperature of 3800K. In K dwarfs, CaH signals are generally 
stronger, and the TiO signature is most prominent in M dwarfs. 

Conclusions. Modeling synthetic polarization signals from starspots for a range of G-K-M dwarfs leads to differences in the 
prominence of various molecular signatures in different wavelength regions, which helps to efficiently select targets and exposure 
times for observations. 
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1. Introduction 


Magnetic fields play a crucial role in the activity phe¬ 
nomena on the Sun and on cool stars. Solar activity fea¬ 
tures have been studied in detail to reveal the structure 
of the underlying magnetic fields, which was used to con¬ 
strain the solar dynamo theory. The detailed physics of 
stellar magnetic interactions, however, remains vague be¬ 
cause we have insufficient direct measurements of mag¬ 
netic fields on the surfaces of cool stars. A comprehensive 
overview of starspot phenomena on different types of cool 
stars, observational tools, and diagnostic techniques for 
studying starspots and their propertie s, such as magnetic 
fields, is given by iBerdvug-i^ ( 2005ll . Stellar magnetic 
fields have been measured i n various stellar objects with 


different techniques (e.g., IValenti fc Johns-Krulll 12001 


Donati fc LandstreetH2009l : iR.einersI 2012 1. Such 


measure¬ 
ments can be achieved directly by the use of polarimet- 
ric techniques. The detection of magnetic fields on cool 
stars via polarimetry, however, is impeded on one hand 
by limited observational tools, that is, there are hardly 


any large telescopes with high-resolution spectropolarime- 
ters. On the other hand, the net circular polarization sig¬ 
nal from regions of mixed polarity fields of the unresolved 
stellar disk cancels out, which only allows detections of 
large-scale magnetic fields. The main part of our knowl¬ 
edge about magnetic fields on cool stars and in starspots is 
based on the analysis of Zeeman-broadened spectral lines. 
Zeeman signatures in atomic lines that are due to mag¬ 
netic fields in starspots are diluted by contributions from 
the rest of the star. To detect such magnetic fields, the 
use of molecular over atomic lines is favorable because the 
latter receive a strong contribution from outside the spot 
umbra, whereas some molecular lines can only be formed 
in cool starspots if the effective temperature of the stellar 
photosphere is high enough. The choice of the adequate 
molecule therefore represents the key to measuring mag¬ 
netic fields unambiguously in unresolved spots wherein the 
molecules are formed. Here, we investigate the four most 
prominent molecular bands of TiO, CaH, MgH, and FeH. 
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Table 1. General comparison of the obtained signals showing the advantages of different molecules for various spectral 
type investigations. Calculations were made with various combinations of underlying models for each spectral type. 
We list the minimum residual intensity (I), Stokes V, and Stokes Q signals from the different molecular lines for a 
model star with vsini=10 km/s and 10% of the longitudinal or transverse magnetic field of 3 kG. Tabulated maxima 
exclude excursions due to atomic lines that overlap the relevant molecular band. 


spectral type 

Tphot/Tspot 

Minimum residual intensity 

Maximum Stokes V/Ic 

[%] 

Maximum Stokes Q/Ic 

[%] 



TiO 

MgH 

CaH 

FeH 

TiO 

MgH 

CaH 

FeH 

TiO 

MgH 

CaH 

FeH 

M 

3500K/2800K 

0.702 

0.671 

0.661 

0.688 

0.707 

0.758 

0.730 

0.704 

0.410 

0.153 

0.414 

0.166 


3500K/3000K 

0.705 

0.667 

0.660 

0.692 

0.988 

1.016 

1.078 

0.829 

0.529 

0.218 

0.599 

0.199 


3800K/3000K 

0.800 

0.690 

0.771 

0.723 

0.659 

0.602 

0.712 

0.623 

0.352 

0.125 

0.395 

0.150 


4000K/3000K 

0.866 

0.701 

0.830 

0.741 

0.498 

0.404 

0.535 

0.513 

0.266 

0.084 

0.297 

0.123 


4000K/3500K 

0.860 

0.692 

0.829 

0.746 

1.086 

1.058 

1.436 

0.877 

0.491 

0.226 

0.616 

0.222 

K 

4200K/3500K 

0.919 

0.715 

0.884 

0.774 

0.833 

0.724 

1.096 

0.729 

0.376 

0.155 

0.470 

0.184 


4500K/3500K 

0.970 

0.806 

0.957 

0.871 

0.600 

0.461 

0.784 

0.576 

0.270 

0.100 

0.336 

0.145 


4500K/3800K 

0.969 

0.800 

0.954 

0.871 

0.693 

0.790 

0.750 

0.790 

0.211 

0.170 

0.395 

0.204 


4800K/3800K 

0.981 

0.912 

0.978 

0.946 

0.534 

0.560 

0.575 

0.651 

0.162 

0.121 

0.303 

0.168 


4800K/4000K 

0.984 

0.906 

0.978 

0.945 

0.461 

0.816 

0.464 

0.798 

0.090 

0.170 

0.268 

0.221 


5000K/3800K 

0.984 

0.954 

0.983 

0.969 

0.459 

0.458 

0.492 

0.580 

0.139 

0.100 

0.260 

0.149 


5000K/4000K 

0.987 

0.948 

0.983 

0.968 

0.397 

0.670 

0.398 

0.712 

0.077 

0.140 

0.231 

0.197 


5000K/4200K 

0.992 

0.941 

0.985 

0.967 

0.248 

0.912 

0.283 

0.804 

0.040 

0.176 

0.168 

0.237 

G 

5500K/3800K 

0.989 

0.987 

0.989 

0.987 

0.334 

0.302 

0.356 

0.456 

0.101 

0.066 

0.189 

0.117 


5500K/4200K 

0.994 

0.978 

0.990 

0.985 

0.182 

0.608 

0.206 

0.634 

0.030 

0.118 

0.123 

0.187 


5500K/4500K 

0.999 

0.977 

0.996 

0.989 

0.035 

0.580 

0.085 

0.471 

0.010 

0.182 

0.041 

0.149 


5800K/3800K 

0.990 

0.992 

0.991 

0.989 

0.285 

0.245 

0.302 

0.403 

0.086 

0.054 

0.161 

0.104 


5800K/4200K 

0.995 

0.984 

0.991 

0.987 

0.156 

0.494 

0.175 

0.562 

0.026 

0.096 

0.105 

0.166 


5800K/4500K 

0.999 

0.983 

0.997 

0.992 

0.030 

0.473 

0.073 

0.419 

0.009 

0.148 

0.034 

0.133 


6000K/3800K 

0.991 

0.993 

0.991 

0.990 

0.258 

0.215 

0.273 

0.374 

0.077 

0.048 

0.145 

0.096 


6000K/4200K 

0.995 

0.986 

0.992 

0.989 

0.141 

0.435 

0.159 

0.522 

0.023 

0.085 

0.095 

0.154 


6000K/4500K 

0.999 

0.986 

0.997 

0.992 

0.028 

0.418 

0.067 

0.389 

0.008 

0.130 

0.031 

0.123 


6000K/4800K 

1.000 

0.992 

0.999 

0.997 

0.003 

0.208 

0.017 

0.174 

0.002 

0.081 

0.008 

0.057 


In contrast to sunspots, our knowledge of starspots and 
particularly of their magnetic properties is very scarce. 
However, over the past years, observational tools and di¬ 
agnostic techniques for studying starspots have continued 
to improve. Long-term measurements (light-curve model¬ 
ing) investigate the changes and periodicity of large-scale 
fields in stellar atmospheres. Mapping the starspot dis¬ 
tribution on the stellar surface was made possible with 
the Doppler imaging technique, which uses high-resolution 
spectral line profiles of rapidly rotating stars. To reveal 
the magnetic field distribution in addition to the tem¬ 
perature and abundance structure on the stellar surface, 
the magnetic Zeeman Doppler imaging method, based on 
the analysis of high-res olution spect r opolarimetric data , 
is the appropriate tool ( Semell [l989t iMorin et al.ll2013h . 
Recently, exoplanet transits, (spectro-) interferometry, mi- 
crolensing, astereoseismology, etc. enlarged the diagnostic 
repertoire that is available for studying starspots. 

The importance of starspots, in addition to studying 
the underlying magnetic field and its evolution, lies in 
their significance for the ongoing exoplanetary research. 
A planet transiting a starspot can reveal parameters of 
the exoplanetary sys t em su ch as its obliquity, as shown in 
Sanchis-Oieda et al. ( 2013h for a number of Kepler planet 
hosts. The presence of starspots outside the path of the 
transit is also relevant because it affects the depth of the 
transit or secondary eclipse from which the size and the 


albedo of the exoplanet can be deduced. To assess this 
effect, parallel observations of the transit or secondary 
eclipse curves as well as of the stellar activity are required. 
For the latter, the properties of the starspot - such as mag¬ 
netic field, temperature, and size - need to be monitored, 
and their effects on the observed signal need to be deter¬ 
mined. 


The technique applied throughout this work is the use 
of molecular lines as direct diagnostics of cool spots on the 
surfaces of active stars. Given a high photospheric tem¬ 
perature, particular molecular lines can only be formed in 
cool starspots, which allows detecting spots irrespective 
of their distribution, and furthermore, on slowly rotating 
stars. 


Vogt ( I 979 I I reported on the first detection of molec¬ 


ular lines from starspots when he found TiO and VO 
bands on a K2 star, indicating the existence of a spot 
with a spectral ty pe of an M6 star, where the mo lecu- 
lar li nes origina t ed. Huenemoerder fc RamsevI ( 1987^ and 
later Neff et aP ( 199,^ an'd O’Neal et al. 1199611 developed 
a technique for determining spot temperatures and filling 
factors from molecular lines. The method fits an observed 
spectrum with the sum of two components containing ef¬ 
fective temperatures of the spot and the unspotted pho¬ 
tosphere weighted by spot filling factor and continuum 
surface flux ratio. Extending this basic approach to the 
polarimetric domain, we model the spectrum of a spotted 
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star using G, K, and M type stars and apply various mag¬ 
netic field strengths to estimate the expected (circular and 
linear) polarization signal. Considering the new generation 
of high-resolution spectropolarimeters, it is important to 
know which molecule is the most magnetically sensitive 
and produces a d e tectabl e signal in circular polarization. 
Berdvugina et al. ( 2n06bll for the first time detected cir¬ 
cular polarization of 0.5-1% in TiO, CaH, and FeH lines 
on three active M dwarfs. 

Detection of Zeeman-broadened or even Zeeman-split 
lines led to detections and estimates of magnetic fields on 
G, K, and M main-sequence stars. As discussed in the com - 
prehensive overview on starspots bv iBerdviigina ( 2005 1. 
approximate values of starspot filling factors, tempera¬ 
tures, and magnetic field strengths are known. In particu¬ 
lar, the models applied in this paper are based on observed 
combinations of spot and photosphere temperatures. 

For stellar surfaces, high filling fa ctors of up to 50% o r 


even more have been determined by lO’Neal et al.l (|l996fl . 


( 200 


(e.g., iNeff et al.Nl995l : lO’Neal et al.Nl99ffl . IValenti et al.l 
(1199811 explored the extent to which the spectral synthesis 
of TiO lines could reproduce the optical spectrum of an in¬ 
active M dwarf. For sunspots, the first s pectropolarimetric 


models and observations were shown in iBerdvugina et al 


This unusual starspot size might be a composition of 
multip le sunspot-sized spots ( Solankil 12002 1. Berdvugin^ 
^ discussed the reasons for the contradiction of such 


( 200 nll . IBerdvugina et al. ( 200.311 concluded that the TiO 
7 ( 0 , 0 )i ?3 system band head at 7054 A is one of the best 
molecular diagnostics of the magnetic field in (sunspot 
umbra and) starspots. They expected a Stokes V signal 
of 0.3% and were able to nicely repro duce TiO line po¬ 
larizat ion observed on three M dwarfs ( Berdvugina et aP 
l2006bll . The triplet states of this system can be described 
sufficiently well by Hand’s case (a), especially for transi - 
tions between lower rotational levels. lO’Neal et all (l2004[l 
presented measurements of starspot temperature and fill¬ 
ing factor on five highly active stars, using absorption 
bands of TiO. They fitted TiO bands using spectra of 
inactive G and K stars to represent the unspotted pho¬ 
tospheres of the active stars and spectra of M stars to 
represent the spots. 

For M gH (and o ther diatomic molecules), iKronid 


large spot filling factors with Doppler imaging results. 

In contrast to the Sun, where large spot areas result 
in a higher luminosity due to a corresponding increase 
in plage area, the photometric and spectroscopic vari- 
ability of active st ars is dominated by the spot umbra 
( Radick et al.lll99f)ll . Low spot filling factors in G dwarfs 
were found. A representative sample of starspot tempera¬ 
tures on very active stars shows a clear tendency for spots 
to show a higher contrast to the photosphere than in hot¬ 
ter stars: the temperature difference between spots and 
the photosphere decreases from about 2000 K in GO stars 
to 200 K in M4 st ars, which holds for dwarfs and giants 
( Berdvuginall2005ll . 

Gooler d warfs have stronger magnetic fields that cover 
larger areas ( Johns-Krull fc Valenti I99d) . To detect mag¬ 
netic fields in the starspot umbra, it is necessary to keep 
track of the spectral lines that are weak outside the spot 
and strong in the umbra; the goal of this work is to reveal 
such molecular lines. 

In this analysis we investigate the temperature range 
in which the molecules MgH, TiO, GaH, and FeH can serve 
as indicators for magnetic fields on very active cool stars, 
and we present synthetic Stokes profiles for the modeled 
spectral type. This study is meant to serve as basis for 
modeling and planning spectropolarimetric observations. 

2. Molecular bands 

TiO lines are prominent in spectra of sunspots and cool 
stars. TiO was first ast ronomically detected in spectra of 
M dwarfs ( Fowleiiri904ll . The first detection of TiO molec¬ 
ular bands fro m starspots was rep o rted by IVogtJ (I 1979 II 
for a K2 star. R.amsev fc Nationa ( 198flll suggested us¬ 
ing TiO bands to measure starspot properties. Further 
observations of TiO bands in spectra of active stars al¬ 
lowed the spot area and temperature to be measured 


(1928) and iHilll ( 1929ll theoretically predicted the molec¬ 
ular Zeeman splitting. Laborato ry measur e ments of the 
Zeeman effect were reported in Grawfortl (Il934h . MgH 
spectral lines were described in lLabordel ( 1961 ). where a 
comp arison of laborato ry wavelengths with solar ones was 
given. Bell et ahl ( 19851) us ed MgH lines to determine the 
surface gravity of Arcturus. IBerdvugina fc Savanov (1992) 
modeled the MgH ra dical as a diagnos t ic of the surface 


gravity of K giants. Berdvugina et al. ( 200Clll presented 


the first spectropolarimetric measurements and modeling 
of MgH lines in sunspots. The A^H state of the MgH A^H- 
(0,0) transition is an intermediate Hund case (a-b) 
with spin-orbit constants of F = 5.7, meaning that this 
system is relatively close to a pure Hund case (b). With 
the Paschen-Back effect taken properly into account, the 
MgH A^H-X^E^ system at 5200 A represents a sensitive 


2OOOI 


tool for stellar magnetic studies (IBerdvugina et al 

I2OO5II . 

Because it is one of the most important astrophysi- 
cal molecules, CaH bands c an also be used as a luminos¬ 
ity indicator of cool stars ( Oehmai] 19^; IMou U 119761: 
Mould fc WallisI 1977 : Barbuv et al.f 


1977: Barbuv et al 


199311 . CaH bands 


represent a valuable tool for studying brown dwarfs since 
the absorption therein is an important opacity source. The 
CaH A^n-X^E+ band system is observed in the wave¬ 
length region at 6600-7600A. The ground state X^E+ 
is well described by a pure Hund case (b). The excited 
state A^H is intermediate between Hund cases (a) and (b). 
Like for MgH, the Pasch en-Back effect must be considered 


(|Berdvugina et al.1 1200311 . The first polarimetric measure¬ 


ments and modeling of CaH St okes profiles was carried 
out bv IBerdvugina et ah ( 2006all . They confirmed the di- 
ag nostic value of CaH i n cool astrophysical sources. 

IValenti et al. ( 200 ill clearly showed the usefulness of 
the FeH lines as stellar magnetic diagnostics with their 
detection of Zeeman broadening of FeH lines in an active 
M dwarf. They also modeled the stellar spectrum using 
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Fig. 1 . Modeling of a G dwarf: Combinations of the photosphere and spot models for a model star with vsini=10 km/s, 
10% magnetic field of 1 kG (solid lines) and 3 kG (dashed lines), Tphot=5800 K, Tspot=4200 K. Stokes V was calculated 
with an inclination angle of the magnetic field 7 = 0° and Stokes Q with 7 = 90°. 


the sunspot spectrum by Wallace et al. ( 1998ll . In these 
sunspot observations, a remarkable magnetic sensitivity 
of many lines of the FeH F'^A-X^A system can be seen. 
It was concluded that this particular molecular system 
would be a powerful tool for diagnosing solar and stel¬ 
lar magnetism once the spin-coupling constants became 


availa ble (e.g., Berdvugina et al. 2000ll . Reiners fc Basri 
(|2006ll measured magnetic fields on ultracool stars using 
FeH and emphasized the fact that its response to mag¬ 
netic fields is easier to detect than other magnetic diag¬ 
nostics. The FeH F'^A-X^A system is produced by tran¬ 
sitions between two electronic states with couplings of the 















































































































N. Afram and S.V. Berdyugina: Molecules as magnetic probes of starspots 


5 


FeTi 


MgH 


Example K dwarf 


FeH 



0.2000 


0.1000 

i ;; 

0.0000 

Of 

lAf 

- 0.1000 


- 0.2000 

■ i , • - 


5200,3 5203.7 

Wavelength [A] 

CaH 

Ti Fe Ti Fe 



Co 


TiO 

Mn 


Cr 


1.002 



1.002 



0.983 



0.983 

- \ \ 

/ V V V V 

0.964 



< 0.964 

- G, / 

- 

0.945 



0.945 

- 

- 

0.926 



0.926 






6936.0 6939.' 

Wovelength [A] 


7055.5 

Wavelength [A] 


Fig. 2. Modeling of a K dwarf: Combinations of the photosphere and spot models for a model star with vsini=10 km/s, 
10% magnetic field of 1 kC (solid lines) and 3 kG (dashed lines), Tphot=5000 K, Tspot=4000 K. Stokes V was calculated 
with an inclination angle of the magnetic field 7 = 0° and Stokes Q with 7 = 90°. 


angular momenta that are intermediate be tween limiting 
Hund cases (a) and (b) (IPhillips et a,l.lll987). Its high mag - 
netic sensitivity was pointed out bv I Wallace et al.l (|1998ll . 
With the help of th e previously missing spin-coupling con¬ 
stants provided by Dulick et al.l (|200,3ll . the perturbation 
calculation of the molecular Zeeman effect for this tran¬ 


sition allowed Afram et al. ( 20071 l2008h to analyze the 
FeH F'^A-X'^A system and its diagnostic capabilities for 
investigating solar and stellar magnetic fields. They pre¬ 
sented fits t o the first spectrop olarimetric measurements 


in sunspots ( Afram et al.l 1200711 . providing Lande factors 


for energy levels and transitions. They also compared syn- 
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Fig. 3. Modeling of an M dwarf: Combinations of the photosphere and spot models for a model star with vsini=10 km/s, 
10% magnetic field of 1 kG (solid lines) and 3 kG (dashed lines), Tphot=3500 K, Tspot=3000 K. Stokes V was calculated 
with an inclination angle of the magnetic field 7 = 0° and Stokes Q with 7 = 90°. 


thetic a nd observed Stokes profiles in different wavelength 2.1. The model 
regions (jAfram et al.11200811 . 

The Stokes parameters I, V, and Q w ere calculated 
with the radiative trans fer code STOPRO ( Solankilll987l : 
Berdyugina et all 200, 3h . which assumes local thermody¬ 


namic equilibrium and solves the radiative transfer equa¬ 
tions upon input of molecular or atomic linelists and an 
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atmosph eric model. For this wo rk we used Phoenix stellar 
models ( Hauschildt et al. 199(tl in local thermodynamic 
equilibrium (LTE) with an effective temperature range of 
2800 K to 6000 K and a surface gravity of 4.5, correspond¬ 
ing to G, K, and M dwarfs. In dwarfs, the diatomic hy¬ 
drides MgH, FeH, and Ca H are more abundant and hence 
stron ger than in giants ( Bell et al. 1985; Barbnv et al.l 


199,111 . However, the magnetic signatures studied in this 


work are independent of the log g value. In contrast, other 
diatomic molecules, such as CN and CO, are stronger in 
giants. The radiative transfer code STOPRO treats the 
spectral line in a magnetized atmosphere accounting for 
the atomic and molecular Zeeman and Paschen-Back ef¬ 
fects. If not declared otherwise, we assumed the star to 
have a circular spot with an area of 10% of the stellar disk 
and a spot comprised of only longitudinal or transverse 
magnetic fields. We obtained the spectrum of the spotted 
star by integrating synthetic polarized spectra over the 
stellar surface taking into account limb darkening, stellar 
rotation, and the direction of the magnetic field. We esti¬ 
mated the strength of the circular and linear polarization 
Stokes V and Q signals for the different molecules MgH, 
TiO, CaH, and FeH. We combined various photosphere 
and spot models according to the c urrent knowledge of 
stellar properties (|Berdvuginal 1200511 . For CaH and TiO, 
a number of atomic lines are included in the calculations, 
as indicated in Figs. [T]131 We took into account the follow¬ 
ing line-broadening effects: intrinsic, thermal (according 
to the local temperature in the model atmosphere), mag¬ 
netic, rotational (10 km/s), microturbulent (1 km/s), and 
instrumental (0.07 A) broadening. 


tries in the polarization signals for MgH and CaH occur 
as a result of the Paschen-Back effect and lead to a net 
polarization signal across the line profiles. In Fig. [I] for 
TiO some strong polarization signals are due to blended 
atomic lines, for example, the Co line near 7054 A. Table 
[T] tabulates the maximum amplitude of the polarization 
signal only from the molecular species, which forms exclu¬ 
sively in the spot. 

Resulting minimum residual intensities, the maximum 
Stokes V and Q amplitudes (in units of the continuum 
intensity) for MgH, FeH, CaH, and TiO are given in Table 
□ for all calculated G dwarfs. 


3.2. K dwarfs 


The effective temperatures of K stars range from 4000 K 
to 5000 K. Spot temperatures Tsnn t on K dwarfs aver¬ 
age 73% of Tphot (IBerdvuginal 1200511 . Spot filling factors 
have values from 20% - 60% and magnetic field filling 
factors range from 10% - 80%, for typical magnetic field 
strengths of 1.5 kC - 2.5 kC, with one exception of 5 kG. 
In Fig. □ results for a single type of K star are presented 
for the molecules MgH, FeH, CaH, and TiO, combining a 
photosphere model with Tphot=5000 K and a spot model 
Tspot=4000 K. For TiO, a broad Stokes Q signal can be 
seen, which is promising to be detected even with low- 
resolution observations. See Table □ for the exact values 
of the residual intensity, Stokes V, and Stokes Q signals 
and for all other calculated K dwarf combinations. 


3. Modeled signals for G-K-M dwarfs 

3.1. G dwarfs 

Type G dwarfs with photosphere temperatures of 5000 K- 
6000 K show an average spot temperature Tgpot of 70% of 
Tphot- There are exceptions, however, such as EK Dra, 
which have a lower contrast between Tphot and Tgpot 
(Tphot = 5900 K, Tspot = 4900 K). Spot filling factors are 
relatively low, ranging from 10% - 25% on active dwarfs 
and from 30% - 35% on active giants. The Sun has a spot¬ 
filling factor of 1%. Magnetic field measurements for the 
Sun show an umbral magnetic field strength of 3 kC. 

In Fig. □ we present synthetic Stokes profiles for the 
molecules we investigated for a G dwarf. For clarity, a 
single type of a G star is presented in the figure, with the 
combination of the photosphere model with Tphot=5800 K 
and a spot model Tspot=4200 K (see Table □for all calcu¬ 
lated combinations). We applied an instrumental broad¬ 
ening of 0.07 A and a rotational broadening corresponding 
to V sin i=10 km/s. 

In each figure, the middle panel (Stokes V) was calcu¬ 
lated with an inclination angle 7 = 0° (longitudinal mag¬ 
netic field) and the lowest panel with an inclination angle 
7 = 90° transverse magnetic field). In the latter case, a 
relatively strong Stokes Q can be seen. Strong asymme- 


3.3. M dwarfs 


The effective temperature for an M dwarf lies between 
2600 K and 3800 K. Spot temperatures tend to be closer 
to Tphot , with Tspot being approximately 86% of Tphot- 
Observed spot filling factors are low with 5% -20%, but 
magnetic field filling factors show high values from 50% - 
85% (see discussion in Sect. 13.411 . The range of observed 
mag netic field strengt hs is rather wide with 2.5 kC - 4.3 
kC ( Berdvugina 200^ . To model an M dwarf (Fig. [3]), we 
combined a photosphere model with Tphot=3500 K and 
a spot model Tspot=3000 K for the molecules MgH, FeH, 
CaH, and TiO. In the TiO spectrum, a broad Stokes Q 
signal is present here as well. In M dwarfs, atomic lines are 
much more blended with molecular lines than in G and K 
dwarfs. We tried to select wavelength regions with a lower 
atomic line contribution than from molecular lines. We 
find that for MgH in M dwarfs, the strength of the Stokes 
V signal for some lines is almost the same, regardless of 
the magnetic field, because of line saturation. In addition, 
the Paschen-Back effect causes nonlinear blending of the 
Zeeman components. Stokes Q is not saturated because it 
has more allowed Zeeman transitions, and, therefore, sat¬ 
uration does not occur yet. We list all minimum residual 
intensities, Stokes V, and Stokes Q signals for all calcu¬ 
lated M dwarf combinations in Tabled 
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Fig. 4. Combination of the photosphere and spot models 
for the MgH A^n-X^I]+ system at 5200 A for a model 
star with vsini=10 km/s, Tphot=5000 K, rspot=4000 K, 
and with (i) left panel: 5%, 10%, and 20% (solid, dashed, 
and dotted lines, respectively) magnetic field of 3 kG (ii) 
right panel: 10% magnetic field of 1 kG, 2 kG, and 3 kG 
(solid, dashed, and dotted lines, respectively) magnetic 
field. Stokes V was calculated with an inclination angle of 
the magnetic field 7 = 0° and Stokes Q with 7 = 90°. 


3.4. Sensitivity of the modeled signal 


To illustrate the sensitivity of the modeled signal to 
changes in spot properties, in Fig. U we have modeled a 
K dwarf for the MgH A^n-X^S+ system at 5200 A with 
three values of the spot size: 5%, 10%, and 20% (in the left 
panel) and with three values of the magnetic field strength: 
1 kG, 2 kG, and 3 kG (right panel). In this regime of 
the Paschen-Back effect, the signal still increases with ap¬ 
parent spot size or magnetic field strength. These spot 
sizes can still be called moderate, giv en that filling fac - 
tors of up to 80% have been reported (|Berdvuginall2005ll . 
In fact, the magnetic field can be complex, and opposite- 
sign-polarity fields cancel out in the polarization signal, so 
that the actual magnetic filling factor can be much higher 
than the assumed 10%. The minimum residual intensity 
reaches 0.96, 0.95, and 0.94, the Stokes V amplitudes (in 
units of the continuum intensity) are 0.32%, 0.67%, and 
1.20%, and the Stokes Q amplitudes (in units of the con¬ 
tinuum intensity) are 0.07%, 0.14%, and 0.22% for spot 
sizes of 5%, 10%, and 20% respectively. For magnetic field 
strengths of 1 kG, 2 kG, and 3 kG, the minimum resid¬ 
ual intensity is 0.95 in all cases, the Stokes V amplitudes 
(in units of the continuum intensity) are 0.34%, 0.53%, 
and 0.67%, and the Stokes Q amplitudes (in units of the 
continuum intensity) are 0.02%, 0.08%, and 0.14%. 


4. Conclusions 


We modeled various scenarios to calculate the expected 
circular and linear polarization signal in MgH, TiO, GaH, 
and FeH lines present in starspots. We used a model as¬ 
suming a spot size of 10 % of the stellar disk comprising 
either longitudinal or transverse magnetic fields and es¬ 
timated the strength of the Stokes V and Q signals in 
molecular lines that could be expected to be observed 
from starspots on G-K-M dwarfs using realistic estimates 
of starspot parameters. For M dwarfs and cooler K dwarfs, 
it should be noted that molecular lines can also be formed 
in the (cool enough) photospheric or nonmagnetic model, 
which is taken into account in this analysis. Galculations 
were made for the model combinations of the spectral 
types G, K, and M, as plotted in the previous section, 
and for a range of G, K, and M dwarfs with different un¬ 
derlying models. In Table [1] the resulting minimum values 
of the residual intensity, V, and Q are listed for all calcu¬ 
lated model combinations and for the following molecular 
lines: the MgH line at 5197.35 A, the (blended) FeH lines 
at 9923.44 A and 9923.35 A, the GaH line at 6937.03 A, 
and the (blended) TiO lines at 7056.43 A. 


An overall comparison of the obtained signal shows 
the advantages of different molecules for various spectral 
type investigations in Fig.[S] In general, the Stokes V and 
Q signals are stronger for cooler spot temperatures be¬ 
cause of the stronger molecular lines. Glear differences 
can be seen in the usefulness of the analyzed molecules 
for the different spectral types and for different contrasts 
of the assumed photosphere and spot temperatures. In 
G dwarfs, the considered MgH and FeH lines show the 
strongest Stokes V and Q signals, while other MgH lines, 
such as the line at 5201.6 A are, in principle, stronger than 
FeH, but are blended by atomic lines. In K dwarfs, MgH 
shows strong Stokes V signals, but the Stokes Q signals 
are generally stronger in GaH. For M dwarfs, GaH shows 
an overall stronger signature an d may be more readily de- 
tected than TiO (as predicted in iBerdvngina et al.l (l2006a^ 
and MgH. TiO itself, however, exhibits strong Stokes V 
and Q signals especially in M dwarfs. The Stokes V sig¬ 
nals for FeH are the strongest in all G dwarfs, with a spot 
temperature of 3800K and in almost all K dwarfs with 
this spot temperature. Atomic lines and other molecular 
lines were not included in the MgH spectral range cal¬ 
culations to solely show Paschen-Back features. Although 
large telescopes with high-resolution spectropolarimeters 
exist, detecting magnetic fields on cool stars remains a dif¬ 
ficult task. To use the available observing time in its full 
capacity, it is important to efficiently select targets and 
exposure times upon prior modeling, as presented in this 
work. 
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Fig. 5. Maximum Stokes V signal (left panel) and Q signal (right panel) from the different molecular lines for different 
spectral types for a model star with vsini=10 km/s and 10% longitudinal or transverse magnetic field of 3 kG. In the 
plot, the applied spot temperatures are given and the photosphere temperatures are implied by the indicated spectral 
types. 
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